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We have studied the antiferromagnetic order - disorder transition occur- 



£> ■ ring at T = in a 2-layer quantum Heisenberg antiferromagnet as the inter- 

im 

plane coupling is increased. Quantum Monte Carlo results for the staggered 

^^ ■ structure factor in combination with finite-size scaling theory give the critical 

ON ■ ratio J c = 2.51 ± 0.02 between the inter-plane and in-plane coupling con- 



stants. The critical behavior is consistent with the 3D classical Heisenberg 



e 

i-rt ■ universality class. Results for the uniform magnetic susceptibility and the 

a: 

correlation length at finite temperature are compared with recent predictions 



for the 2+1-dimensional nonlinear u-model. The susceptibility is found to 
exhibit quantum critical behavior at temperatures significantly higher than 
the correlation length. 



It was recently suggested |T|-§[ that the unusual normal-state magnetic properties of 
the high-T c superconducting cuprates are characteristic of two-dimensional (2D) quantum 
antiferromagnets close to the critical point of a zero-temperature order - disorder transition, 
with the disordered phase having a gap towards spin excitations. It has been argued that the 
physics of such antiferromagnets is described by the nonlinear cr-model in 2+1 dimensions ||. 
Studies of this field theory based upon a 1/iV expansion have resulted in detailed predictions 
for the behavior of near-critical systems |]-f| • In order to test these predictions, it is useful 
to compare them with exact numerical results for some appropriate model. The 2-layer 
Heisenberg antiferromagnet can be tuned through an order - disorder transition by varying 
the coupling between the planes |7||| , and constitutes an ideal system for such comparisons. 
In this Letter, the T = order - disorder transition and the finite-temperature "quantum 
critical" regime of this model are studied using a modification of the Handscomb quantum 
Monte Carlo algorithm PJlQf. Details of this work will be presented elsewhere [[Ll| . 

The model we study is defined by the hamiltonian 

H = J\ 2J Z2 ^<M ' Sa,j + J2 z2 ^U ' ^2,i (1) 

a=l,2 (ij) i 

where (i, j) is a pair of nearest-neighbors on a square lattice, and S a> i is a spin-| operator 
at site i in plane a. With the inter-plane coupling J 2 = 0, the independent planes have 



long-range order at T = fl2|l , and the spectrum is gapless. For a large ratio J = J2/J1, 
there is a tendency for neighboring spins in adjacent planes to form singlets. There is a gap 
for spin-1 excitations and no long-range order. A series expansion calculation by Hida gave 
a critical coupling J c = {J2I J\) c — 2.56 0. A Schwinger boson mean-field calculation by 
Millis and Monien, on the other hand, resulted in J c = 4.48 ||. 

The coupling ratio J is analogous to the coupling g of the 2+1-dimensional nonlinear 
cr-model. In their study of this model, Chakravarty et al. || identified three regimes in 
the T — g plane. For g < g c there is long-range antiferromagnetic order at T = 0. At low 
temperatures, in the so called renormalized classical (RC) regime, the correlation length £ 
diverges as e 27rps,/T , where p s is the spin-stiffness. For g > g c , there is an excitation gap 



and the correlation length is constant in the low-temperature "quantum disordered" (QD) 
regime. For g « g c , £ ~ T _1 in the high-temperature "quantum critical" (QC) regime. 
Exactly at g c , p s vanishes and the QC regime extends down to T = 0, whereas for g ^ g c 
there is a cross-over to either the RC or the QD regime as the temperature becomes low 
enough for the deviation from g c to be sensed. On the lattice, the spins become effectively 
decoupled as T —> oo and there is a high-temperature cross-over from the QC regime to a 
"local moment" (LM) regime. 

The 3D nonlinear cr-model is the appropriate continuum field-theory for the phase tran- 
sition of the 3D classical Heisenberg model. The T = transition of 2D quantum antifer- 
romagnets is therefore expected to belong to the universality class of that model, provided 
that the cr-model description is valid at the critical point ||. 

Chubukov et al. 0j3|] showed that close to criticality, many physical observables depend 
in a universal manner on a few model-dependent parameters. Once these parameters are 
determined, the temperature dependence of e.g. the wave-vector and frequency dependent 
magnetic susceptibility is known for temperatures T < J\. 

Quantum Monte Carlo studies have confirmed that the 2D Heisenberg model has long- 
range order at T = [I2|. The low-temperature behavior is consistent with the predictions 



for the RC regime [|13|,|14|] • It has been argued that this model is close enough to criticality to 
exhibit QC behavior for 0.35 < T/Ji < 0.55 |2|||. However, this regime is narrow, making 
it difficult to verify the predicted behavior. Introducing frustrating interactions reduces the 
long-range order and widens the QC regime. Unfortunately, frustrated quantum models are 
difficult to study numerically, due to "sign problems" which arise in Monte Carlo algorithms 



15 1 . The 2-layer model (1) does not have this problem, and can be tuned through the 
critical point by varying J 2 / ' J\. 

In order to determine the critical ratio J c = {Jij ' J\) c of the 2-layer model, and to inves- 
tigate its T = critical behavior, we have carried out quantum Monte Carlo simulations of 
periodic lattices with 2L 2 spins, with L = 4, 6, 8, 10. In order to obtain essentially ground 
state results we chose an inverse temperature j3 = 48, which for the system sizes studied 



is sufficient for all calculated quantities to have saturated at their T = values. Monte 
Carlo moves necessary to ensure ergodicity in the subspace with zero total magnetization 
[J2a,i Sai = 0] were carried out. We have also investigated the finite temperature properties 
for various values of J near J c . In these finite-temperature simulations, Monte Carlo moves 
changing the total magnetization were carried out. Systems with L up to 24 at T/Ji > 0.3 
were studied [|X6 | . For small systems we have checked simulation result against exact di- 



agonalization data. At higher temperatures our results are in good agreement with series 
expansion results recently obtained by Singh and Sokol |L7|. 



We have calculated the in-plane staggered structure factor for coupled L x L planes 

Si(L) = ^(sU+iShK-i) 1 ^ (2) 

i,l 

and the full two-plane staggered structure factor 

S *( L ) = ± E([^M + / - Wft - SyX-l)^. (3) 

i,l 

In addition we have evaluated the corresponding staggered susceptibilities Xi an d Xii with 

1 ? 

Xi(L) = I5 E/^^i lW (r)^(0))(-l)^ (4) 

i,l Q 

and a similar expression for \i- 

Two possible order parameters of the phase transition are the sublattice magnetizations 
m\ and m,2 of a single plane and the whole system, respectively. These can be defined in 
terms of the structure factors as 

m n (L) = yj3S n (L)/nL 2 . (5) 

For J < J c the asymptotic T = spin-spin correlation functions 

Ci(f)= (Sl l+r Si t )(-iy* +r y (6a) 

<? 2 (r> ([Sl i+r - Sl t+r ][Sl - SlM^r +ry (6b) 

should have the form 



C n (r)=m 2 n + b n r-^\ (7) 

which gives for the sublattice magnetization 

ml(L) = m 2 n (oo) + k n {\jhf-\ (8) 

Exactly at the critical point, we expect that 77 is equal to the 3D Heisenberg exponent 
7] k, 0.03 ||18| . Hence, we have fit our results for w? n to (|]) with this rj. For m\ the Monte 



Carlo results agree well with this form for all L > 4, whereas L > 6 is needed to obtain good 
fits to the results for m\. Fig. 1 shows mf(L) versus 1/L for J = 2.4, 2.5, and 2.6, along with 
least-squares fits of (§) to the J = 2.4 and 2.5 data. At J = 2.5 the extrapolated values of 
mi(oo) and 777,2(00) are both zero within statistical errors, indicating that the critical ratio 
is close to 2.5. 

Define a reduced coupling j = (J — J c )/J c . As j —* from above, the correlation length 
£ diverges as j~ u , and the staggered structure factors and susceptibilities diverge as j~ 7s 
and j~ 7x , respectively. These exponents are related according to 

75 = „(1 - v ) (9a) 

Tx = 1/(2-17). (9b) 



For a quantity A which diverges as j 7A , finite-size scaling |T9[ relates the value Al for a 
finite system to the infinite-size value Aoo according to 

A L (j)=A 00 (j)f[UU)/L}- (10) 

Eqs. (H) and (|T0| ) give for the size-dependence of S n (L) and Xn(^) at the critical point: 

S n (L,j = 0)~L 1 - r > (11a) 

Xn (L,j = 0)~L 2 -". (lib) 

Fig. 2 shows results for hi(S n ) and ln(x„) versus hi(L) at J = 2.5. If Eqs. ( |TTD hold, the 
data should fall onto straight lines with slopes 1 — rj and 2 — rj, respectively. All of the Si 



results agree well with this form, whereas the other quantities agree within statistical errors 
for L > 6. 

In order to test whether the exponent v agrees with its expected 3D Heisenberg value 
v « 0.70 | IBf , one can use the scaling relation ( |I0"D for j > 0. Graphing Al(j)P a versus Lj u 



for various J and L should produce points collapsed onto a single curve. This is indeed the 
case for S n and Xn if J c ~ 2.50. The best over-all results are obtained with J c = 2.51, in good 
agreement with the results for mi. Fig. 3 displays results for Si using v = 0.70, 77 = 0.03, 
and J c = 2.51 for various J and L. Based on the appearance of such graphs with various 
assumptions for J c , and the results for mi displayed in Fig. 1, we estimate the critical 
coupling and its error limits to be J c = 2.51 ± 0.02. This is only slightly lower than Hida's 
series expansion result (J c = 2.56) 0. 

We now discuss some finite-temperature results for systems close to criticality. For the 
single-plane Heisenberg model, QC behavior has been observed for the uniform susceptibility 
at temperatures 0.35 < T/Ji < 0.55 @|§. For the 2-layer model with J « J c , one would 
expect the cross-over from the LM regime to the QC regime to occur at a higher temperature, 
since the thermal fluctuations are reduced by the strong coupling between the planes. 

Chubukov et al. || carried out -^ expansions of the non-linear cr-model and obtained the 
temperature dependence of a number of observables. Their result for the uniform magnetic 
susceptibility is (for N = 3) 

X, = ^n(^l)(^ +0 .79 3 7XT), (12) 

where c is the spin-wave velocity. Hence, at the critical coupling, where p s = 0, the suscepti- 
bility graphed versus the temperature should produce a straight line with intercept zero and 
a slope which depends only on c. Fig. 4 shows numerical results for the q = susceptibility 

Xu = nTjBKi + s^Mj + SSj]) (13) 

for L=10. The size- dependence of Xu is very weak for L > 10 at the temperatures studied. 
For J = 2.5 a least-squares fit to the T < 0.9 results gives an intercept close to zero. The 
slope of the line gives c = 2.39 ± 0.02, in units of Jia/Ti [a is the lattice constant]. 



The inverse correlation length is predicted to be a linear function of the temperature: |3j 

rl . L0791x21ll( ^_i£|. (14) 



In order to extract £, we fit the correlation function Ci(f), Eq. fl6a]), to the form Ci(r) = 
A e -r/€ r -(i+v) with 77 = 0.03. We have taken the effects of the periodic boundaries into 
account analogously |lTJ to the proposal for ID systems in Ref. |[20| . In Fig. 5, L = 10 



and L = 24 results for £ _1 at J = 2.5 and 2.6 are graphed along with the predicted form 
([UP for p s = and c = 2.39. In view of the susceptibility results one would expect (|i~4"D 
to apply for T < 0.9, which is clearly not the case. It appears that a cross-over from the 
LM regime occurs around T = 0.6. Below this temperature the J = 2.6 results are close to 
the predicted form. For J = 2.5 the correlation length grows faster than 1/T in the regime 
studied, suggesting that this coupling is slightly lower than the critical coupling. Larger 
systems at lower temperatures have to be studied in order to determine the T — > behavior 
of £ more precisely. 

In conclusion, we have studied the order - disorder transition of a 2-layer Heisenberg 
antiferromagnet using a quantum Monte Carlo technique |9|{n3]. The critical ratio between 
the inter-plane and in-plane coupling constants was determined to be 2.51 ± 0.02. The 
T = critical behavior is consistent with the transition belonging to the universality class 
of the 3D classical Heisenberg model. At finite temperature we have studied the uniform 
magnetic susceptibility and the correlation length. Close to criticality the susceptibility is a 
linear function of the temperature for T < 0.9Ji, in agreement with predictions 0,[| for the 
2+1-dimensional nonlinear cr-model. The predicted linear behavior of the inverse correlation 
length applies only for T < O.6J1. 

For the single-plane 2D Heisenberg model, the uniform susceptibility has a roughly linear 
temperature dependence for 0.35 < T < 0.55 [f2]|| . However, a regime where £ _1 has the 
T dependence predicted for the quantum critical regime has not been observed fi"3| , |2~T|| . 
The results presented here indicate that the uniform susceptibility exhibits the predicted 
linear behavior at temperatures significantly higher than the correlation length. Hence, it 



appears that the susceptibility exhibits quantum critical behavior well beyond the cross-over 
boundaries defined by the behavior of the correlation length. 
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FIGURES 

FIG. 1. The sublattice magnetization mi versus 1/L for J = 2.4 (solid squares), J = 2.5 
(open squares), and J = 2.6 (solid circles). The dashed and solid curves are least squares 
fits of the form given by Eq. (^ with rj = 0.03 for J = 2.4 and 2.5 respectively. 

FIG. 2. Size dependence of Si (open squares), S2 (solid squares), xi (open circles), and 
X2 (solid circles) at J = 2.5. The solid and dashed lines have slopes 1 — 77 = 0.97, and 
2 — r] = 1.97, respectively. 

FIG. 3. Finite-size scaling of S\ with J c = 2.51 and 3D Heisenberg exponents. Open 
squares are for J = 2.55, solid squares for J = 2.60, open circles for J = 2.70, solid circles 
for J = 2.75, and crosses for J = 2.80. The solid curve is the asymptotic form for the scaling 
function, f(x — ► 0) ~ x 1_?? . 

FIG. 4. The uniform susceptibility versus the temperature for L = 10 at J = 2.4 (open 
squares), J = 2.5 (solid squares), and J = 2.6 (open circles). The line is a least-squares fit 
to the T < 0.9, J = 2.5 data. 

FIG. 5. The inverse correlation length for L = 10 (open circles) and L = 24 (solid circles) 
versus the temperature for J = 2.5 and J = 2.6. The solid lines are of the predicted form, 
Eq. (HD with p s = and c = 2.39. 
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